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ABSTRACT 

We study a sample of 39 massive early-type lens galaxies at redshift z < 0.3 to determine the 
slope of the average dark-matter density profile in the innermost regions. We keep the strong 
lensing and stellar population synthesis modeling as simple as possible to measure the galaxy total 
and luminous masses. By rescaling the values of the Einstein radius and dark-matter projected 
mass with the values of the luminous effective radius and mass, we combine all the data of 
the galaxies in the sample. We find that between 0.3 and 0.9 times the value of the effective 
radius the average logarithmic slope of the dark-matter projected density profile is — 1.0 ± 0.2 
(i.e., approximately isothermal) or —0.7 ±0.5 (i.e., shallower than isothermal), if, respectively, a 
constant Chabrier or heavier, Salpeter-like stellar IMF is adopted. These results provide positive 
evidence of the influence of the baryonic component on the contraction of the galaxy dark-matter 
halos, compared to the predictions of dark matter-only cosmological simulations, and open a new 
way to test models of structure formation and evolution within the standard ACDM cosmological 
scenario. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: structure — dark matter 



1. Introduction 

Understanding the original formation and the 
subsequent evolution of the galaxies we observe 
today remains one of the major open questions 
in modern astrophysics. Within the past decade, 
there has been incredible progress in the realiza- 
tion of high-resolution numerical simulations that 
are now starting to reproduce in detail the phys- 
ical properties of the galaxies obse r ved at red- 



shift z 







2007: Tissera et al 



Meza et all 120031: iNaab et all 
20101) . iV-body cosmological 



simulations have predicted that in an expanding 
universe cold dark-matter particles collapse into 
gravitationally bound, self-similar halos with a di 



verging inner density p rofile (e.g., iNavarro et al 



19961 : iMoore et allll999h . The values of the three- 
dimensional logarithmic slope 7 = dlnp/dlnr of 
the collapsed dark-matter halos have been found 



to be approximately equal to —1 and —3, respec- 
tively, in the innermost and outermost regions. It 
is in these halos that the stars of the observed 
galaxies were assembled. Recent hydrodynami- 
cal simulations have shown that several mecha- 
nisms associated to baryonic physics affect the 
stellar mass assembly of a galaxy (e.g., dissipa- 
tionless accretion of stars originally formed far 
from a galaxy center and dissipational gas flow- 
ing towards the inner regions of a galaxy, later 
transformed into stars). The complex interplay 
between the luminous and dark components can 
alter significantly the dark-matter distribution in 
the center of a halo, making it steeper or shal- 
lower, depending on the role played by the diffcr- 
ent physical processes (f or more details, see e.g. 



Lackner fc Ostrikerll2010D . 
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In the last two decades strong gravitational 
lensing combined with stellar dynamics and/or 
stellar population synthesis models has been ex- 
tremely successful in measuri ng; the amount and 
distribution of dark matter ( Grillo et al.l 12008c, 
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20101: iGrillo fc Christen"seiill201ll:lBarnabe fc Koopmansl 



20071 iBarnabe et all 120091 : lAuger et al.1 1200.,, 
2010al: iFadelv et al.l |2010|), the presence of dark 



DPI 

9, 



matt e r substructure (e.g ., Vegetti fe Koopmansl 



matt er halos (e.g. 
20101 : 



2009(; IVegetti et al"1l201C), and the sizes of dark 



Suvu fc Halkolal2010HRichard et al 
Donnarumma et al.l 1201 lF ) in early-type 
galaxies beyond the local Universe. The com- 
bination of these mass diagnostics has also al- 
lowed to find alternative ways to address some 
interesting astrophysical and cosmological top- 



initial mass function (IMF) (e.g., Grillo et al. 


2008a. 2009: Gril 


o & Gobad 2010: llreu et al. 


2010: Auger et al. 


2010bt Spiniello et al.l 201 ll 


Sonncnfcld ct al. 2011) and of the values of 



the co s mological paramet e rs (e. g. . IGr illo ct al 



Schwab et al 



200814 IParaficz fc Hjorthl 120091 : 
20101 : ISuvu et al.ll20irl) . 

The Sloan Lens ACS (SLACS) survey has been 
crucial for the identification of a statistically sig- 
nificant sample of strong gravita tional lensing sys- 
tems. Disparate studies (e.g.. | Treu et al.l 1200 " 



Bolt on et al.l 120061 : IGrillo et al.ll2009t lAuger et al. 
2009) have shown that the SLACS lens galaxies 
are a representative sample of the parent sam- 
ple of massive early-type galaxies observed in 
the Sloan Digital Sky Survey (SDSS). By mod- 
eling the strong gravitational features detected 
in these lensing systems, it has been possible 
to obtain accurate and precise total mass esti- 
mates pro jected w i thin t h e corresponding E i nstein 
radii llTreu et al.l 120061 iKoopmans et all 120061: 
Bolton et all 120081 lAuger et al l l2009h . Here, we 
exploit the fact that the Einstein radius of a lens- 
ing system is not a length scale intrinsic to the lens 
(since it depends also on the redshift of the source) 
to study the average inner dark-matter density 
distribution of a specific lens sample. We do this 
by combining the lens aperture total and lumi- 
no us mass measurem ents. The pioneering work 
of iRusin et al. ( 2003 ) prefigures to some extent 
the general method and results presented here. In 
this previous analysis a self-similar mass model 
for early-type galaxies was constrained by using 
aperture mass-radius relations from 22 gravita- 
tional lenses. The total mass distribution of the 
lens galaxies was described in terms of a two- 
component (luminous and dark matter) model 
parametrized by (1) a present-day normalization 



value of the i?-band stellar mass-to-light ratio, (2) 
the dependence of a galaxy S-band stellar mass- 
to-light ratio on its luminosity, (3) the projected 
dark over total mass fraction within two effective 
radii and (4) the three-dimensional logarithmic 
density slope of the dark-matter profile. 

This Letter is organized as follows. In Sect. 
2, we introduce the sample of massive early-type 
lens galaxies. In Sect. 3, we describe the method 
and hypotheses used to determine the inner slope 
of the average galaxy dark-matter density pro- 
file. In Sect. 4, we illustrate the main results of 
this analysis. In Sect. 5, we compare our results 
with those of previous studies and anticipate fu- 
ture prospects. In Sec. 6, we draw conclusions. In 
the following, we assume Ho = 70 km s _1 Mpc -1 , 
fl m = 0.3, and = 0.7. 

2. The sample 

In this work, we concentrate on 39 mas- 
sive early-type lens galaxies discovered in the 
SLACS survey and studied in several papers (e.g., 
Bolton et aLll2008l: iGriilo et al.ll2009l: lAuger et ail 
2010a|). In detail, we conservatively consider only 
those galaxies that satisfy the photometric and 
spectrosc opic selection criteria of the sample an- 
alyzed in iGrillol (j2010l ) (i.e., values of the SDSS 
fracDeV morphological index larger than 0.95 in 
the r, i, and z bands; SDSS spectroscopic redshifts 
z sp between 0.05 and 0.33; SDSS aperture stellar 
velocity dispersions between 150 and 400 km s — 1 ; 
total luminous masses between 10 10 5 ans 10 12 
Mq). These galaxies have both accurate total Mt 
and luminous Ml mass estim ates obtained from , 
respectively, strong lensing ( Auger et al.l l2009h 
and spectral energy distribution (SED) fitting 
models (from the public galaxy catalogs provided 
by the MPA/JHU collaboration^). The physical 
properties of the galaxies in the sample are sum- 
marized in Table [TJ This is a specific sample of 
early-type galaxies with large values of central 
stellar velocity dispersion a® (more details on the 
measure ments of the physical quantities can be 
found in lGrilldl2010h . Therefore, the results of 
the analysis performed in this letter should not 
be simplistically generalized to early-type galaxies 
with different physical properties until verified by 
larger samples. 



x http: / /www. mpa-garching.mpg.de/SDSS/ 
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Table 1: Physical properties of the early- type lens galaxies of the sample. 



•^sp Re 


^Ein 


^0 


M L 


M T (< R Ein ) 


(kpc) 


(kpc) 


(km s- 1 ) 


(10 10 M ) 


(10 10 M ) 


0.06-0.32 3.2-16 


1.3-7.0 


200-320 


7.4-56 


3.9-47 



Notes - Ranges of values of the spectroscopic redshift z sp , effective radius i? e , Einstein radius Rmn, central 
stellar velocity dispersion ao, luminous mass Ml (assuming a constant Chabrier stellar IMF), and total 
mass projected within the Einstein radius Mx(< i?Em)- 

References - SDSS and MPA/JHU public catalogs: IXueer et al.1 (l2009h : iGrilld (<2010h . 



3. The method 

For the lenses in the sample, we measure here 
the values of the dark-matter mass density pro- 
jected within the Einstein radius and study, in a 
statistical way, their dependence on the projected 
distance from the lens centers. 

In practice, we proceed as follows. For each 
lens galaxy, we define an adimensional radius A 
as the ratio between the Einstein radius i?Ein and 
the effective radius R P : 



A 



Rev 



(i) 



The Einstein radius of a lens galaxy depends on 
its total mass distribution, but also on the redshift 
of the lensed source. Thus, i?Ein is not a funda- 
mental property of a galaxy and we use A instead 
to quantify the distance from the center of a lens. 
The latter is a scale-free distance that is obtained 
by normalizing the value of the Einstein radius to 
the typical scale of a galaxy luminous mass distri- 
bution (i.e., R e ). 

Similarly, we estimate the value of the dark- 
matter mass projected inside the cylinder with ra- 
dius equal to the Einstein radius as the difference 
between the values of the total Mt{< Rkiu) and 
luminous Ml(< R E - ln ) masses and rescale the re- 
sult to the total amount of luminous mass Ml of 
each galaxy. We notice that the measurements of 
the projected masses within R E [ n are robust and 
almost model-independent for the total ones and 
only scaled according to the fraction of total light 
of a de Vaucouleurs profile for the luminous ones 
fsee lGrillo et~aDl2009l lGrilldl2010h . Then, we de- 
fine an adimensional dark-matter projected mass 
density ^ as the ratio between the adimensional 



value of the dark-matter projected mass and the 
area of the disk with radius equal to the value of 
A: 



* : = 



M T (< R Ein ) - M L (< R Ehl ) 1 



Mr 



ttA 2 



(2) 



In this way, both A and $ are referred to the lumi- 
nous properties of the galaxies in the sample and 
can thus be properly compared!! 

Next, we measure the value of the Kendall 
rank correlati on coefficient g (for its definition, see 
ISalkindl 2006) and check if the values of A and \1/ 
are correlated at a statistically significant level. In 
the case of a significant correlation, we perform a 
Markov chain Monte Carlo study on the galaxy 
sample to characterize the joint probability distri- 
bution function of the values of the two coefficients 
a and j5 that are used to fit a power-law relation 
to our set of data: 



* = a x (Af . 



(3) 



We apply this method to the sample described 
m Sect. 2, starting from different hypoth eses. 



First, we assume a constant Chabrier (2003) (la- 
beled as Ch) stellar IMF to estimate the lumi- 
nous mass values of all the galaxies in the sam- 
ple. Then, we rescale the ga laxy lum i nous mass 
values to a constant heavier, ISalpeterl (Il955h -like 
(labeled as Sa) stellar IMF by simply multiplying 
the Chabrier luminous mass values by a constant 
factor equal to 1.7. Next, we consider the case of a 



2 In pa ssing, we notice that differently from iRusin et al.l 
(2003) the luminous mass values of the lens galaxies are 
measured here from the multi-band photometric and spec- 
troscopic observables and are not scaled according to the 
galaxy B-band luminosity values. 
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non-universal stellar IMF and mimic a variation, 
moving from a lighter to a heavier IMF (labeled as 
Ch — > Sa), depending on the values of the galaxy 
central stellar velocity dispersion. This is moti- 
vated by the facts that stellar velocity dispersion is 
currently considered the most significant parame- 
ter related t o the stellar populat ion properties of a 
galaxy (e.g., Graves et al. 2009) and that a stellar 
IMF variation with stellar velocity dispersion has 
been tentatively detected bv lTreu et al. ( 2010l ). In 
detail, following the previous indications, we use 
a toy model in which we multiply the Chabrier 
luminous mass values with a factor that increases 
linearly from 1.0 to 1.5 as the value of ao changes 
from 200 to 320 km s" 1 . 

We conclude by remarking that the correlation 
of the errors on A and ^ is not significant and 
therefore will not affect our results on the steep- 
ness of the average dark-matter density profile. 
Although obtained from the same sets of observa- 
tional quantities, the uncertainties on A are very 
small (the median relative error is smaller than 
4%) and mainly related to the quality of the pho- 
tometric measurements, while the uncertainties on 
^ are considerably large (the median relative er- 
ror is approximately 40%) and primarily driven by 
the degeneracies that are inherent in the popula- 
tion synthesis modeling. 

4. Results 

We summarize in Table [2] the values of the 
Kendall rank correlation coefficient g between A 
and VP and remark that using the three different 
hypotheses mentioned above about the stellar IMF 
of the sample galaxies always results in an anti- 
correlation of the values of A and f at a statis- 
tical significance level higher than 97%. In the 
same table, we also show the best-fitting (mini- 
mum chi-square) /?best and the 68% CL interval 
$68% cl values of the average inner slope of the 
dark-matter projected mass density. These num- 
bers are obtained from Monte Carlo chains with 
5 x 10 5 points for each of the three cases. The data 
set and the best-fitting power-law for the case of 
a constant Chabrier stellar IMF are illustrated in 
Fig. [1] and the marginalized probability distribu- 
tion functions of f3 for the three cases are plotted 
in Fig. [2] 

From Tableland Figs. [T]and [2j we notice that 



^ 1.0 




Fig. 1. — The adimensional values of the dark- 
matter projected mass density within the Einstein 
radius, ^P, and Einstein radius, A. The points, 
with their 1 a error bars, are obtained by using 
the values of the total luminous mass and effective 
radius of the galaxies as dimensional scales and as- 
suming a constant Chabrier stellar IMF. The best- 
fitting power-law is shown in gray. 

assuming a constant Chabrier stellar IMF leads to 
an average dark-matter density profile that consid- 
ered in terms of a three-dimensional spherical pro- 
file decreases in the inner regions approximately as 
1/r 2 , i.e. like an isothermal profile (as usually re- 
ferred to in lensing studies). A constant heavier 
stellar IMF results in a broader probability distri- 
bution function for /3, centered on a slightly larger 
value. This result can be qualitatively explained 
in the following way. If we keep the values of the 
total mass fixed and increase those of the lumi- 
nous mass, we obtain values of ^ that are on av- 
erage smaller and decrease less steeply with in- 
creasing values of A than in the previous case (see 
Fig. [TJ. This translates into a dark-matter den- 
sity profile that is shallower than an isothermal 
one in the center. On the contrary, the proposed 
variation in the stellar IMF provides a steeper pro- 
file of the dark-matter component in the inner re- 
gions. This result can also be understood looking 
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Table 2: Correlations and power-law fits of A and vF 



e (A,*) 


Pbest 


A}8% 


CL 


Ch -0.57 (< 0.01) 


-1.04 


[-1.26, 


-0.78] 


Sa -0.24 (< 0.03) 


-0.77 


[-1.14, 


-0.15] 


Ch -> Sa -0.52 (< 0.01) 


-1.28 


[-1.54, 


-0.93] 



Notes - Values of the Kendall rank correlation coefficient g between A and W (in parentheses, the probability 
that an equal number of measurements of two uncorrelated variables would give values of the coefficient 
higher than the measured ones), and of the best-fitting /?bcst and 68% CL interval /?68% cl of the inner 
slope of the average dark-matter projected mass density. 



at Fig. Q] As expected, the values of ao are pos- 
itively correlated with those of A. This follows 
from the fact that more massive galaxies yield, on 
average, larger Einstein radii. Therefore, varying 
the stellar IMF from a Chabrier to a Salpeter-like, 
the points in Fig. Q] with small values of A have 
approximately the same values of $ (because of 
the unchanged Chabrier stellar IMF), while those 
with large values of A have now larger luminous 
mass values (because of the changed, heavier stel- 
lar IMF), hence, in general smaller values of ^F 
The net effect is an increase in the value of the 
slope (3. 

We notice that the dark-matter universal profile 
obtained from dark matter-on ly cosmological sim- 
ulations ( Navarro et al.lll996h is characterized by 
values of j3 of approximately —0.1 and —0.2 within 
0.1% and 1% the value of the typical dark-matter 
length scale (r s ), respectively. 

5. Discussion 



We compare here our results with those of sev- 
eral other studies on early-type galaxies and indi- 
cate a possible way to extend this work. 

Based on a sample of 16 massive Coma galax- 
ies, with physical properties very similar to those 
of the galaxies in our sample. iThomas et al. (2011) 
find th at if the stellar IMF is universal and lKroupa 
(|200lh -like. i.e. very similar to a Chabrier IMF, 
then the galaxy dark-matter density profiles are 
smooth and on average close to isothermal out to 
several tens of kiloparsecs (see Fig. 6 in the cited 
paper). This conclusion follows from joint dynam- 
ical (Schwarzschild's orbit superposition) and stel- 
lar population models that exploit accurate photo- 
metric and spectroscopic data. Recalling that the 



luminous mass estimates obtained by assuming a 
Chabrier or a Kroupa stellar IMF are only slightly 
different, the findings of our study on the value of 
(3 in the 'Ch' case are consistent with those of the 
analysis performed in the Coma cluster. 

iNapolitano et al.l (|201dh consider a sample of 
335 local early-type galaxies and estimate their lu- 
minous and total masses from, respectively, pho- 
tometric SED fitting and dynamical Jeans model- 
ing. They also adopt a Kroupa stellar IMF and 
conclude that the average three-dimensional loga- 
rithmic slope 7 of the dark-matter density profile 
at R e ranges between —2.1 and —1.7. In the sim- 
plified case of a spherical power-law density profile, 
the values of j3 and 7 are related in projection in 
the following way: j3 = 7+ 1 (if 7 is different from 
— 1). Our estimates of (3 in the 'Ch' case are there- 



fore co nsistent with the results of lNapolitano et al 
(|2010h . 



In t he gravitational lensing study bv lRusin et al 
(|2003l ). detailed in Sect. 1, the authors come to 
the conclusion that it is not possible to measure 
precisely the slope of the dark-matter component 
because of the significant degeneracies between 
the parameters. Despite that, models with a 
dark-matter density profile that is approximately 
isothermal are generally preferred to models with 
shallower dark-matter density distributions (e.g., 
with 7 = —1). Our results confirm these last 
findings. 

In the last few years, some new observational 
constraints have been obtained on the stellar IMF 
of massive early-type galaxies. For this specific 
class of galaxies, if the IMF is constant , a Salpeter- 



like IMF is favored by the dat a (e.g. , IGrillo et al 



2008aL 120091 : IGrillo fc Gobatl 120101: iTreu et al 
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1.0 -0.5 0.0 
P 



0.5 



Fig. 2. — Probability distribution functions of the 
average logarithmic inner slope, /3, of the dark- 
matter projected mass density. The different his- 
tograms refer to the different assumptions on the 
stellar IMF discussed in the text. As a refer- 
ence, the arrow close to the x-axis shows the result 
that a three-dimensional spherical density distri- 
bution decreasing as 1/r 2 (i.e., an isothermal pro- 
file) would give. Larger and smaller values of (3 
correspond, respectively, to shallower and steeper 
profiles with respect to an isothermal one. 



20101 ; lAueer et all l2010tj : ISpiniello et "all l201ll) . 

In the case of a constant Salpeter IMF, our 
results indicate a dark-matter density profile 
that is shallower than an isothermal one in 
the central regions Inte restingly, the study of 
Treu fc Koopmansl (|2004 ) on the average inner 
power-law slope of the dark-matter halos of 5 
early-type lens galaxies at z sp « 0.5 — 1.0 pro- 
vides also indication of profiles shallower than 
isothermal ones. The slope values are robustly de- 
termined by combining gravitational lensing and 
stellar dynamics, with and without priors on the 
lens stellar mass-to-light ratios from the Funda- 
mental Plane. More recently, ISonnenfeld et al 



(|201ll ) have also performed a two-component lens- 
ing and dynamics analysis to decompose the total 
mass distribution of the double Einstein ring grav- 



itational lens in terms of a bulge of stars and a 
dark-matter halo. They find that a Salpeter IMF 
is preferred to a Chabrier IMF for the stellar com- 
ponent and that the value of the three-dimensional 
logarithmic inner slope 7 of the dark-matter halo 
is — 1.7 ± 0.2. Therefore, our findings in the 'Sa' 
case are in general good agreement with these 
lensing and dynamics analyses. 

The results of the two combined strong lens- 
ing, s tellar dynamics, a nd st ellar population stud- 
ies bv lTreu et all (|2010l) and lCardone et all (|201ll ) 
on samples of more than 50 SLACS lens galaxies 
agree on finding that a constant heavy (Salpeter- 
like) stellar IMF requires a shallower dark-matter 
density profile than a constant light (Chabrier- 
like) stellar IMF. Furthermore, based on differ- 
ent s amples of SLACS l e nses and m ass diagnos- 



tics, 



Jiang fc Kochanek ( 20071 ) and lAueer etlrt" 



(|2010bt ) conclude that adiabatically compressed 
models of the galaxy dark-matter halos are fa- 
vored. These findings are also in qualitative agree- 
ment with our results. 

The next natural step towards a clearer pic- 
ture of the internal structure of massive early- 
type galaxies will be the extension of the SLACS 
sample to the lens galaxies selected from the 
BOSS (Baryon Oscill ation Spectroscopic Survey; 
lEisenstein et alj 2011 ) Emission- Line Lens Survey 
(BELLS; iBrownstein et"ai]|2012j) . The lens galax- 
ies for which strong lensing and stellar population 
models will be available at the end of this new 
survey will allow to enlarge significantly the lens 
sample (in nearly the same luminous mass range) 
and to explore the average density profile of the 
galaxy dark-matter halos on a radial range (i.e., 
A) that is approximately twice as large as done 
here. 

6. Conclusions 

We have combined strong gravitational lensing 
and stellar population synthesis models in a ho- 
mogeneous sample of massive early-type galaxies 
to measure the logarithmic inner slope of the aver- 
age dark-matter density profile. We have obtained 
clear indication of the contraction of the halos 
when compared to the results of dark matter-only 
cosmological simulations. This is in line with the 
recent findings of high-resolution hydrodynami- 
cal simulations which include radiative cooling 
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and feedback proce sses (e.g.. lAbadi et al.l l20ld 
Tissera et all l20ld : iDuffv et all l2010h . These 
studies show that the contraction of a halo does 
depend not only on the amount and distribution 
of the baryonic mass condensed at the halo centre, 
but also on the details of the halo assembly his- 
tory. Future theoretical and observational efforts 
towards a better understanding of the inner dark- 
matter structure and the stellar initial mass func- 
tion of galaxies will therefore be crucial to explore 
different cosmological models and to investigate 
the nature of dark matter and its interaction with 
baryons. 

C. G. is grateful to Marco Lombardi, Giuseppe 
Bertin, Matteo Barnabe, and Simona Vegetti for 
interesting discussions. This research was sup- 
ported by the DFG cluster of excellence "Origin 
and Structure of the Universe" . 
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